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Abstract: Clustering algorithms build jets though the iterative application of single par- 
ticle and pairwise metrics. This leads to phase space constraints that are extremely com- 
PLi plicated beyond the lowest orders in perturbation theory, and in practice they must be 

^ , implemented numerically. This complication presents a significant barrier to gaining an 

^ analytic understanding of the perturbative structure of jet cross sections. We present a 

I— ' novel framework to express the jet algorithm's phase space constraints as a function of clus- 

pg tered groups of particles, which are the possible outcomes of the algorithm. This approach 

K*" highlights the analytic properties of jet observables, rather than the explicit constraints 

^vq on individual final state momenta, which can be unwieldy at higher orders. We derive the 

0^ form of the n-particle phase space constraints for a jet algorithm with any measurement. 

We provide an expression for the measurement that makes clustering effects manifest and 
f-^ relates them to constraints from clustering at lower orders. The utility of this framework 

CN is demonstrated by using it to understand clustering effects for a large class of jet shape 

• . observables in the soft/collinear limit. We apply this framework to isolate divergences 

. ^H and analyze the logarithmic structure of the Abelian terms in the soft function, providing 

rS the all-orders form of these terms and showing that corrections from clustering start at 

Cd next-to-leading logarithmic order in the exponent of the cross section. 
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1. Introduction 

Jet algorithms play an essential role in high energy collision experiments, organizing the 
hadronic structure of an event to allow for a meaningful interpretation of the short distance 
interaction. The LHC has increased interest in understanding and utilizing the structure 
of jets [1, 2], and these tools often make use of clustering algorithms such as anti-kx, 
Cambridge/Aachen (C/A), and kx [3, 4, 5, 6, 7]. Understanding the structure of jets 
requires understanding their perturbative description. 

Clustering algorithms build jets through an iterative recombination procedure that 
sequentially identifies the pair of particles that are closest according to a distance measure 
and then clusters them, i.e. combines their momenta, doing so until a stopping criterion 



is reached. Each algorithm differs mainly in the choices of the distance measure and the 
stopping criterion used in this procedure. Implementing phase space constraints from 
clustering for final states with more than two or three particles is very complicated, and 
can present an impasse to analytic perturbative calculations, which in many cases must be 
performed numerically. This makes it difficult to understand the effects of clustering on 
the perturbative series in jet cross sections. 

1.1 The Complication of Measuring Properties of Jets 

When measuring a set of observables, {O}, clustering can strongly affect their values by 
pulling particles into and out of the various regions where {0} is measured. When calcu- 
lating these effects on the differential cross section, it is useful to separate the calculation 
into two parts - the squared matrix element, A, and the measurement function, A4, which 
implements the measurement of {O} on the final state. The measurement function depends 
on the final state momenta at a given point in phase space, ^, and the observables being 
measured. It is integrated against A when evaluating the contribution to the cross section: 



da 

d{d} 



I d<^A{<^)M{{0},'^). (1.1) 



We will be concerned with understanding the general structure of the measurement function 
for a wide range of observables for jets defined by a clustering algorithm. This requires 
translating the iterative recombination procedure of clustering algorithms into all-orders 
properties of the measurement function - it is a nontrivial task. 

The measurement function is composed of two main parts: a function of momenta, 
fri{k}), associated with the measured observable Or G {0} and a restriction, TZr{{k},^), 
on the region of phase space over which we measure it. Here {k} is a subset of all momenta 
for a given phase space point, $, that TZr requires to be in region r. The measurement 
function can be expressed schematically as 

M{{0},^) =ll6(Or-Y^fr{{k})nr{{k};^)^ , (1.2) 

{k} 

where the product is over the regions of phase space, r, which contribute to the cross 
section. For example, the constraints that contribute to a dijet cross section are illustrated 
in Fig. 1 and correspond to the two jets, TZj-^ and TZj^, and an out-of-jet region, T^outi 
where a different function of momenta, fr{{k}), may be measured in each region. It is 
often simpler to work with the Fourier transform of the measurement function with respect 
to the observables: 



/oo 
dO, e-"''-^'-7W({0},$) 
-oo 



JJexp I -iXr Y,fr{{k})nr{{ky,^) ) . (1.3) 

{k} 




Figure 1: Phase space regions TZr over which observables Or are measured, contributing to a dijet 
measurement function. 

An illustrative example is the measurement of the invariant mass of a jet, m'j; it 
contributes a factor to Ai of 

M{m^J,^)=6im^J-J2k^■kj^ZJ{{ki,kJ};^)\ , (1.4) 

where Tlj[{ki,kj};^) gives the phase space constraints for particles i,j to be in the jet 
J. Expressed explicitly, TZj is given by a product of theta functions of momenta in $ 
that constrain fcj and kj to be in the jet. In general, whether ki and kj are in the jet 
depends on the clustering procedure of the jet algorithm, meaning TZj depends on all final 
state momenta <5. While the form of Eq. (1.2) is compact, it hides correlations between 
particles, since TZj depends on momenta not in the jet. When the measurement function 
is used in fixed order calculations, the hidden phase space constraints from clustering must 
be made more explicit. However, the form of Eq. (1.2) makes the perturbative structure 
of the jet cross section difficult to understand analytically without explicitly carrying out 
higher order calculations. 

The JADE algorithm [8] provides an example of why it is important to make manifest 
the effects of clustering that are hidden in TZr in Eq- (1-2)- Although JADE is a well- 
defined infrared (IR) safe algorithm that combines particles based on an invariant mass 
metric, clustering effects at higher orders change the structure of the leading log series 
(a" In ") and spoil resummation at this order. This fact is obscured by expressing the 
action of the algorithm solely through the iterative application of its metric. It required an 
analytic calculation at two loops [9] to understand the perturbative structure of jet cross 
sections defined with JADE. More recently it has been shown that clustering effects for 
the C/A and kx algorithms change the structure of the perturbative series at 0(a^ln ) by 
explicit calculation of the Abelian terms, likely spoiling resummation at this order. Such 
effects are absent for anti-kx, providing a theoretical preference for this algorithm. 

These examples illustrate the importance of having an analytic understanding of the 
perturbative structure of jet cross sections. However higher order analytic calculations 
are difficult and can be made prohibitively complex by clustering phase space constraints. 
Since the structure of the measurement function, Ai, determines clustering effects, we 
focus on expressing A^ in a form that makes these effects manifest. This will provide 



insight into the perturbative structure of jet cross sections without having to carry out 
exphcit calculations. 

1.2 Reformulating the Constraints from Clustering 

Running a clustering jet algorithm divides up phase space into mutually exclusive regions 
that cover all of phase space. For an A^-jet observable theses regions may correspond to each 
jet, TZjj^, . . . , T^jjv) and an out-of-jet region T^out- The jet algorithm does this by separating 
final state particles into groups that characterize the outcomes of the algorithm. For 
example, for three particles with momentum ki, k2, and ^3, the possible outcomes, lCi{^), 
of an algorithm are 

• ki, k2, and k^ all cluster together: /C3 

• ki and k2 cluster, but k^ does not (+ permutations): /Ci2,/Ci3, and /C23 

• None of the particles cluster: /Ci 

Each outcome JCi{^) represents a set of phase space constraints requiring the particles 
to cluster into a set of groups. Summing over all possible outcomes, i, of a clustering 
algorithm at a given number of final state particles gives a unitarity relation: 

l = Y^]Ce{<^). (1.5) 

Since the ICi makes the phase space constraints for a given outcome of the clustering 
algorithm explicit, TZr is simplified. TZr{{k}] ^) constrains {k} to be in region r where the 
observable Or is measured, but it no longer must implement the constraints from clustering, 
which appear in the /C^. This means TZr depends only on the momenta in the region r: 

7^,({A:};$) ^ 7^,({A:}) . (1.6) 

To write the measurement function, we simply weight Eq. (1.5) by Eq. (1.2): 

A^({o},«i>)=j;/cK<&)n'^(a-E^(W)^'-(W))- (1-7) 

e r {k} 

While it is simple to enumerate the outcomes of clustering, each /C^ contains a com- 
plicated set of phase space constraints; as such, in order to explore this structure we need 
to know where the IR divergences arise, since they can contribute large logarithms in the 
perturbative series of the cross section. This depends on both the squared matrix element, 
A, and the measurement, Ai. Although the structure of ^ is complicated at higher orders, 
we will see that examining the structure of Ai alone can provide insight. The differential 
cross section is sensitive to the IR regions of phase space determined by TZr in Eq. (1.3). 

In this paper we provide a framework to use unitarity relations of the type Eq. (1.5) 
to express the outcome of the clustering algorithm and isolate the most divergent regions 
of phase space, which contribute the most logarithms in the perturbative expansion. As 
we shall see, this is done by using unitarity iteratively to rewrite exclusive constraints as 
inclusive ones. For example, the case where ki and /c2 cluster but ^3 does not contains a 




Figure 2: Phase space for particle k^ for fixed fci,fc2- In the region /C12, fca would be clustered 
with fci and k2. The exclusive constraint ^(fcs ^ /C12) is shown in (a) by the orange region and the 
inclusive constraint 1 — ^(fcs € /C12) is shown in (b) by the orange region. 

phase space constraint excluding k^ from some region of phase space; it cannot be close 
enough to ki or k2 to fall in one of the other outcomes. This is shown in Fig. 2(a), where 
for a given fci and ^2 there is a region /C12 that k^ can not occupy in order to avoid 
being clustered^. Divergences in phase space arise as fcg" , fc-j' — t- along the axes and the 
associated logarithms will depend on ki and k2 in a nontrivial way through the boundary 
of /Ci2- A more useful description of clustering would rewrite this in terms of an inclusive 
phase space constraint shown in Fig. 2(b), where the exclusive and inclusive ways to write 
this constraint would be 



exclusive : 6{hi /C12) , 

inclusive : 1 - e{k3 G /C12) . (T8) 

The function 6(x) is the usual step function, giving 1 if x > and otherwise. The 
inclusive description expresses the divergences in k^ in a way that is independent of ki and 
/c2 and subtracts a finite contribution from the region /Ci2. As a result, integrating over 
the momenta of particle k^ with the inclusive constraints simplifies the divergent structure 
[10, 11]. Rewriting the exclusive constraint as an inclusive one uses a (simple) unitarity 
relation describing the phase space constraints between particles. One can iteratively apply 
unitarity relations of this kind to rewrite the outcome of clustering only in terms of inclusive 
phase space constraints. 

1.3 An Application: Soft Clustering and Jet Shapes 

Energetic particles have the ability to change the number, energy, and direction of jets. 
Often, though, energetic particles are confined to small regions of phase space and only 



^We decompose the four-momentum as 



k^ — n ■ k h n ■ k h fcf , 

2 2 ^' 



and define k^ = n ■ k and k 
with h^ — 1. 



k. Here n and n are light-cone vectors, n*^ = (1, h) and n^ = (1, — n) 



cluster within jets. This occurs for jet observables where the measurement of an observable 
O constrains the collinear radiation to be tightly collimated within its own jet, compared 
to the jet size; e.g. measuring the jet mass, m,j, when mj/Ej <C R forces this to be true. 
This is satisfied for a range of jet shapes and substructure observables [12, 13, 14, 15, 16] 
in the soft /collinear regime. In this case the boundaries of the jets are determined by 
clustering amongst soft particles. 

In this regime the function of momenta that we measure, /r({^}) in Eq. (1.3), simplifies 
to be a function of momentum of a single particle fr{ki) for a large class of observables. 
For the example of jet mass given in Eq. (1.4), the expression greatly simplifies to give 

Y,h-kjnj[{h,kj}) -^ ^pj-kinj{ki), (1.9) 

where pj is the sum of the collinear momentum and is conserved for each jet direction in 
the soft/collinear regime. The measurement function still contains correlations between 
particles through the outcomes of the algorithm, given by ICf 



X({x},cI>)=J^/C,($)nexp(-iE^-^ 
e r \ i 



{ki)'Rr{ki) . (1.10) 



Soft clustering for the anti-kx algorithm is particularly simple: jets in the soft/collinear 
regime are circular with jet size R around the collinear direction up to power suppressed 
corrections [7, 13, 17]. In this case the sum over outcomes in Eq. (1.10) simplifies to a 
single outcome: none of the particles cluster, and the constraint is trivial, /C($) = 1. This 
removes all correlations between different particles in the final state from the measurement 
function. As a result the measurement simplifies greatly: 

M^).^{{x},<^) =\{M^^\{x},h) . (1.11) 

i 

where 

A^W ({x}, k)=e^^(-i^xr fr{k) 7^^(A;)) . (1.12) 

r 

It is well known that the measurement function also takes this simple form for a large class 
of event shapes in the soft/collinear limit, such as thrust, heavy jet mass, angularities, and 
N-jettiness [18, 19, 20, 21], which makes such shapes theoretically appealing. Notice that 
jet rates do not satisfy the form in Eq. (1.11) except for fixed cone algorithms, which have 
known split/merge problems [22, 23]. 

Algorithms such as C/A or kx modify the phase space regions from the simple cones 
of the anti-kx algorithm, due to clustering amongst soft gluons. This occurs when a soft 
particle is pulled in or out of the jet, and we refer to this as 

boundary clustering : clustering of soft particles across the jet boundary of radius R 
around the jet direction. 

Example configurations for 2 and 3 particles are shown in Fig. 3. Since boundary clustering 





Figure 3: Configurations of two (left) and tliree (rigiit) particle boundary clustering that changes 
the jet boundary. The final state particles shown in black cluster, and the red particle is the result. 
In the two particle case the particle outside the jet is pulled in by clustering, and in the three 
particle case the particles in the jet are pulled out. 



modifies the region over which the measurement is made and is an effect that intrinsically 
depends on the momenta of all particles that are boundary clustered, the phase space con- 
straints for the outcome, /C^($), of the C/A or kx measurement functions will depend on 
multiple momenta in the final state. This introduces correlations into the measurement 
function, which in this case is of the form Eq. (1.10). A useful way to organize the cor- 
relations introduced hy ICi{^) is through a correction to the measurement function of the 
form 

n 

M^^\{x},'^)=llM^'\{x},h)+Mi±{{x},'^), (1.13) 

where M^^' is given in Eq. (1.12) and n is the number of final state particles. This is 
obtained from Eq. (1.10) by a trivial application of the unitarity relation ICi{^) = 1 — /C^ (<&) 
for each outcome i, where this relation defines ICi{^). 

Clustering for algorithms other than anti-kx affect the logarithmic terms arising from 
soft divergences. In the soft and collinear regimes, boundary clustering can lead to large 
logarithms in the cross section, starting at 0{a'l). This was shown in [24, 25] and further 
explored in [26, 27], and the original authors termed the effect clustering logs. In [27] it 
was shown that clustering logs have the same properties as non-global logarithms (NGLs) 
[28, 29] in the soft/collinear limit and that they contribute at least at next-to-leading 
logarithm^ (NLL), but could in principle contribute at LL. It was also argued that new 
clustering effects likely arise at each order in a^ at NLL and beyond, spoiling resummation. 

Using unitarity relations of the form Eq. (1.5), we can express the measurement func- 
tion to isolate the divergences from boundary clustering. These corrections can be pa- 
rameterized in terms of the number of particles involved in clustering. In particular we 



^We count in the exponent of tlie distribution, where N*LL terms are of order a" In' 



n-fc + l 



show 

where AA^''^' is the contribution where k particles cluster with at least one other particle 
and first arises at 0{a^). This result applies not only to boundary clustering, but clustering 
in general, away from soft-collinear regime. While clustering effects can a priori contribute 
at LL, the form of the measurement function in Eq. (1.14) allows us to show that clustering 
logarithms arise at NLL and to determine to all orders the form of the Abelian soft function 
and show that exponentiation of these terms does not occur. 

1.4 Outline 

In this paper, we construct a framework to describe the phase space constraints from 
clustering and the effect on observables. The formalism we use is very general, applying 
to a wide class of clustering algorithms and measurements. It can be applied to a specific 
type of clustering, such as boundary clustering, or to clustering in general. The formalism 
to describe clustering is presented in two main parts. 

First, in Sec. 2, we define a set of functions that schematically describe the phase 
space constraints from clustering. They are functions of clustered groups of particles, 
which are the possible outcomes of the jet algorithm. The physical interpretation of these 
functions allows us to construct unitarity relations that transform exclusive constraints 
into inclusive ones. The result is a prescription to systematically determine the phase 
space constraints from clustering in terms of inclusive constraints. Our goal is not to write 
clustering constraints explicitly in terms of the momenta of individual particles, but instead 
to express the constraints in a form that highlights the analytic properties of the outcome 
of a generic clustering algorithm. We do not deal directly with the individual clustering 
steps of the algorithm, but in most cases our results can be applied to describe each step. 

Second, in Sec. 3 we apply this framework to the measurement function. For any 
measurement function of the form Eq. (1.7) where the observables depend on the momenta 
of a single clustered group, we use the results of Sec. 2 to find the corrections due to 
clustering. We show that these have the form Eq. (1.14). This form can be exploited to 
determine the organization of logarithms arising from clustering. We also comment on the 
infrared safety of clustering algorithms in Sec. 4; we can use this framework to see precisely 
how clustering effects among soft particles are IR safe. 

In Sec. 5, we apply the formalism of the measurement function to study the logarithmic 
structure of the soft function arising from clustering. We focus on the Abelian terms, 
where the matrix elements factorize and it is straightforward to translate the measurement 
function into contributions to the cross section. We use the form of the measurement 
function along with constraints from renormalization to determine the all-orders structure 
of the soft function, allowing us to show that clustering NGLs contribute at most at NLL 
at all orders in a^. Finally, in Sec. 6, we summarize our results. 



2. Phase Space Constraints from Clustering Algorithms 

In this section and the next, "clustering" can refer to general clustering or to a specific kind, 
such as boundary clustering. If this framework is used to describe boundary clustering, 
which is relevant for observables that constrain the collinear radiation in the jet to be 
collimated compared to the jet size (e.g. mj/Ej <C R), then the phase space constraints 
will ignore clustering that does not change the jet boundary. We will also treat the particles 
as indistinguishable; this is not necessary, but simplifies the discussion. 

As discussed in the introduction, the deterministic nature of jet algorithms allows the 
action of the jet algorithm to be characterized in terms of a set of mutually exclusive 
outcomes, /C^, as in Eq. (1.5). When we run the jet algorithm, sets of particles will cluster, 
and we use these sets to characterize the outcomes of the algorithm. For example, if there 
are 5 particles in the final state labeled qi,...,q5, then a possible outcome from a jet 
algorithm J' is 

J ■{qi,q2,q3,q4,q5} — > {{qi,q2},{q3,q4},{q5}}- (2.1) 

In this case particles 1 and 2 have been clustered, as have particles 3 and 4, and particle 5 
does not cluster. 

This example shows that the action of a jet algorithm is to partition the final state 
particles into groups. We refer to partitions like the one above as labeled partitions. There 
is a one-to-one association between labeled partitions and outcomes of the algorithm. How- 
ever, if the particles are indistinguishable, there are many labeled partitions that describe 
the same basic outcome. In the case of Eq. (2.1), the important characteristic is that 
particles were clustered in groups of 2, 2, and 1; a distinct clustering that is equivalent to 
this outcome is {{gi, 94}, {'?2) ^sji {^s}}- Therefore, we can simplify the labeled partitions 
and describe them in terms of the number of particles clustered in each group. We will 
refer to these simply as partitions. The labeled partition in Eq. (2.1) maps to the partition 
{2,2, 1}. For 5 particles, the possible partitions are 

{5}, {4,1}, {3,2}, {3,1,1}, {2,2,1}, {2,1,1,1}, {1,1,1,1,1}. (2.2) 

In general, we will denote partitions by P = {pi, . . . ,Pk}, where the pi are natural 
numbers. The number of labeled partitions that are equivalent to the partition V is 

where n is the total number of particles, and di is the number of elements of V equal to i. 
The first factor in N{V) is the usual multinomial coefficient, and the second factor accounts 
for the degeneracy in dividing particles among equal sized elements in the partition. 

2.1 Unitarity Relations and Phase Space Constraints 

We will now discuss how we can more specifically determine the phase space constraints 
ICi, defined in Sec. 1, for each outcome, £, of the jet algorithm (which we associate with a 



partition). In words, a partition represents the phase space constraints: 

{pi, . . . ,pk} : pi particles cluster, p2 particles cluster, . . ., Pk particles cluster, 

and these groups of particles do not cluster with each other. (2.4) 

Therefore, we must define two basic functions that act on the elements of a partition: one 
requiring a set of particles to cluster and one requiring different groups of particles not to 
cluster. These can be defined in terms of objects in the partition: 

c{p) : p particles are clustered, (2-5) 

and 

s({pi, . . . ,pk}) : the clusters oi pi, . . . ,pk particles do not cluster. (2.6) 

For completeness, we define s{{p}) = 1 and c(l) = 1 (i.e., no constraint). The function 
s{{pi, . . . ,Pk}) defines an exclusive constraint; for example, it excludes the set of particles 
Pfc from some region of phase space so that they do not combine with the sets of particles 
Pi, . . . ,Pk-i- For a given partition V, the phase space constraints we assign to it are 



Hcip) 
pev 



s{V) . (2.7) 



If we sum over all possible outcomes of the jet algorithm, we obtain a unitarity relation for 
phase space: 



l = Y.N{V) \{c{p) 



V 



pev 



s{V) , (2.^ 



where the sum is over partitions. Each term in the sum corresponds to an outcome fCi of 
the algorithm, so that Eq. (2.8) is a more specific version of Eq. (1.5). 

Some basic examples are useful. For two particles, the partitions are {1,1} and {2}, 
and the unitarity relation is 

l = c(2) + s({l,l}). (2.9) 

This relation defines s({l, 1}). For three particles, the partitions are {1, 1, 1}, {1,2}, and 
{3}, and the unitarity relation is 

1 = c(3) + 3c(2)s({l, 2}) + s({l, 1, 1}) . (2.10) 

Consider the function s({l, 2}). Let us label the single particle qi and the pair of particles 
that cluster (72 and (73. The function s({l,2}) allows gi to go anywhere except for near 
52 and 53. However, we can use a unitarity constraint to rewrite s({l,2}) in terms of an 
inclusive constraint: 

s({l, 2}) = {qi goes anywhere) — {qi would cluster with q2, qs, or q2 + qs) . 

= l-m({l,2}). (2.11) 
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Figure 4: On the left, a configuration of 3 particles for the C/A algorithm that contributes to 
c(2) m{{l, 2}) but not to c(3). In this case particles 92 and q^ cluster first to form the merged pair 
92 +937 shown in red. Particle qi is too far away to cluster with 92 + 93) as it lies outside the 
solid circle of radius R centered at 92 + (Zs- Therefore the 3 particles do not cluster, and so do not 
contribute to c(3). On the right, a schematic picture of the region of phase space for particle qi 
defined by m{{l, 2}), which is the complement of s({l, 2}). 



This relation defines the merging function m({l,2}), which requires gi to be near the 
clustered pair. However, this function does not require that all three particles would 
cluster when running the jet algorithm: 



c(2)m({l,2})/c(3) 



(2.12) 



For example, for fixed momenta of 52 and ^3, the single particle qi could be close enough 
to q2 but sufficiently far from q^ so that all three particles would not cluster. In Fig. 4, we 
show a configuration that is part of c{2)m{{l, 2}) but not part of c(3), as well as a schematic 
region of phase space for m({l, 2}). Only in the case of two particles does m{{l, 1}) = c(2). 
In general, merging requires that at least one particle from each group would cluster with 
at least one particle in another group, such that all the groups are connected. The general 
merging function is 



■m{{pi, . . . ,pk}) '■ the clusters of pi, . . . ,pk particles merge. 



(2.13) 



Our general goal will be to rewrite the exclusive s{V) functions in terms of the inclusive 
clustering and merging functions that require particles to be grouped together. For each 
function s{V), we can define a unitarity relation to rewrite s{V) in terms of clustering and 
merging functions. Note that rewriting s(V) does not change the outcome of the jet algo- 
rithm, but instead rewrites the phase space constraints from clustering in the measurement 
function in a way that makes the perturbative structure of the cross section clearer. 



11 



We introduce some diagrammatic notation which allows us to express these unitarity 
relations in a simple way. These diagrams are helpful in thinking about higher order particle 
configurations and are useful in determining unitarity relations. Instead of clustering in 
3-space, it is sufficient to think of clustering along a line. Each particle will be depicted 
as a dot, and if particles are clustered then we circle them. If the particles do not cluster, 
we put a vertical bar between them. If particles are not on the same line, then their phase 
space constraints do not depend on each other. Finally, when we require groups of particles 
to merge, we circle those groups and place a vertical bar between the groups. In Fig. 5, we 
give the diagrams for the 2 and 3 particle cases. 



n = 2 



■^ 



-•+• 



=(2) s({l,l}) 



n = 3 



where 



= cg> + @ |* + • hh 

c(3) 3c(2)s({l,2}) s({l,l,l}) 

^\^ = :^ - ^^ 

c(2)s({l,2}) c(2) c(2)m({l,2}) 



Figure 5: Clustering diagrams for the cases n = 2 and n = 3. Particles are represented by 
dots, and clustered particles are circled. Particles that do not cluster have a vertical bar between 
them. Merged groups of clustered particles are circled, with a vertical bar between each group. 
The different terms in the unitarity relations are shown. We also show how to express the s{V) 
functions requiring particles to not cluster in terms of c{p) and m{'P) functions that require the 
particles to be clustered/merged. 

2.2 Clustering Maps 

When we rewrite the phase space constraints given by s{V) in Eq. (2.8), we can repre- 
sent the result as a map between partitions. Consider the unitarity relation that defines 
s({l,l,I}) in Fig. 5: 

s({l,l,l}) = l-c(3)-3c(2)s({l,2}) 

= l-c(3) -3c(2) + 3c(2)m({l,2}). (2.14) 
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Each term represents a different phase space constraint on the three particles, and corre- 
spondingly a map to a different partition. 



1 


{1,1,1} - 


■> {1,1,1} 


-3c(2) 


{1,1,1} - 


^ {1,2} 


3c(2)m({l,2}) 


{1,1,1} - 


^ {1,2} 


-c(3) 


{1,1,1} - 


^ {3} 



{3} <"'' 



Each of these maps represent a coarsening of the partition {1,1,1}. A partition Vi is 
coarser than a partition V2 if "Pi is made by combining elements of 7^2 ("^2 is /iner than 
Vi). These maps can be used to determine all ways to rewrite a given 3(1^). For a given 
partition V, we write all maps to coarser partitions (or itself), and each step in the map is 
associated with a factor that either merges or clusters particles. For instance, in the map 
{1, 1, 1} — )• {1, 2} — )• {3}, the associated factors are 

map: {1,1.1} "» {1.2) -* {3} 

contribution: — 3c(2) — m(|l,2}) 

In the first step, we cluster two single particles, and there are 3 ways to choose a pair 
of particles. In the second step, we merge the single particle with the clustered pair. 
The minus signs arise because we are using a unitarity relation, as in Eq. (2.14). Note 
that whenever a map involves only grouping single particles, e.g. {1,1,1} — )■ {1,2}, a 
clustering factor is generated, and whenever it involves a group with multiple particles, 
e.g. {1,2} — ;■ {3}, a merging factor is generated. 

We can think of a map, (p, as a sequence that successively coarsens the initial partition: 
H'Pi^'Pn) = {Pi, . . . , Vn} , with Vi+i coarser than Vi . (2.17) 

Each step along the sequence is associated with a clustering or merging factor. We define 
the contribution from each step in the map to not include the functions s{Vi) appearing in 
the full phase space constraints. To formalize this, we will define a few special partitions: 






the base partition {1, . . . , 1} . (2-18) 

the partition that gives the outcome of the algorithm 
the final partition in a map (p . 



In Eq. (2.8), we obtain a factor of s{V^) that we want to rewrite in terms of inclusive c 
and m functions. The sum over all possible maps cl){V^,V-'), including summing over all 
possible V^ , gives the unitarity relation that lets us rewrite s{V^). By definition, V-' must 
be coarser than or equal to V^. 

As an example, for 4 particles we can use these maps to determine the unitarity relation 
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for s({l, 1, 1, 1}) = ■s(7^'^). These maps and the associated factors are 



</): V^ 



{1,1,1,1} 

{1,1,2} 
{1,3} 
{1,3} 
{2,2} 
{2,2} 

{4} 

{4} 

{4} 

{4} 

{4} 

{4} 



Summing over ah maps gives the unitarity relation 



■po 






— ;> 


•po 






— > 


•po 






— > 


-pO 


— >• 




{1,1,2} 


•pO 






— > 


pO 


— >• 




{1,1,2} 


•pO 






— > 


pO 


— )• 




{1,3} 


pO 


— > 




{2,2} 


pO 


— > 




{1,1,2} 


pO 


— )■ 


{1, 


1,2} ^{1,3} 


pO 


— )■ 


{1, 


1,2} ^{2,2} 



contribution 

1 

-6c(2) 
-4c(3) 

[-6c(2)][-2m({l,2})] 

-3c(2)2 

[-6c(2)][-c(2)] 

-c(4) 

[-4c(3)][-m({l,3})] 

[-3c(2)2][-m({2,2})] 

[-6c(2)][-m({l,l,2})] 

-6c(2)][-2m({l,2})][-m({l,3})] 

[-6c(2)][-c(2)][-m({2,2})] 

(2.19) 



s({l, 1,1,1}) = 1 -6c(2) -4c(3) + 12c(2)m({l,2}) -3c(2)2 + 6c(2)2 -c(4) 
+ 4c(3)m({l,3}) + 3c(2)2m({2,2}) + 6c(2)m({l,l,2}) 
- 12 c(2) m({l, 2}) m({l, 3}) - 6 c(2) V({2, 2}) . (2.20) 

We give the clustering diagrams for 4 particles in Fig. 6, which can alternatively be used 
to derive this relation. 

This map structure is very useful once we add in the measurement on the final state. 
There is another important identification to make in these maps. We can divide a map 
p«) _i, p/ jj^^Q ^^Q parts: 



pi 



clustering 



'pc 



merging/clustering 



p/ 



as m 



(j){v^,vf) = M'P'",^') ■ (kc{v',vf) 



(2.21) 



(2.22) 



This will provide a key simplification in the measurement function. The map V^ — )• V^ 
involves only clustering functions. The map V'^ — )■ V^ involves merging or clustering 
functions, but the first step in the map must be merging so that V^ is unambiguously 
defined. V^ can be equal to V^ or 7'-' , depending on the map. To make the role of V^ 
clear, we repeat the 4 particle maps with V^ explicitly identified. 
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n — 4 



c(3) 



^ 



c(2)^ 



c(2) 



c(4) 4c(3)s({l,3}) 3c(2)2s({2,2}) 



+ 



^f 



♦4«- + • • • • 



6c(2)s({l,l,2}) s({l, 1,1,1}) 

c(3)m({l,3}) c(3)s({l,3}) 

+ -©to 



c(2)2m({2,2}) c(2)2s({2,2}) 



c(2)m({l, 1, 2}) 2 c(2)m({l, 2})s({l, 3}) 

c(2)s({l,l,2}) c(2)2s({2,2}) 



Figure 6: Clustering diagrams for the case n = A. The meaning of the diagrams are the same as 
in Fig. 5. These can be used to derive the unitarity relations for 4 particles, such as s({l, 1, 1, 1}). 
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<t> 



qyW 






pc 




-pf 


{1,1,1,1} 






{1,1,1,1} 




{1,1,1,1} 


{1,1,1,1} 






{1,1,2} 




{1,1,2} 


{1,1,1,1} 






{1,3} 




{1,3} 


{1,1,1,1} 






{1,1,2} 




{1,3} 


{1,1,1,1} 






{2,2} 




{2,2} 


{1,1,1,1} 


— ^ 


{1,1,2} 


^ {2,2} 




{2,2} 


{1,1,1,1} 






{4} 




{4} 


{1,1,1,1} 






{1,3} 




{4} 


{1,1,1,1} 






{2,2} 




{4} 


{1,1,1,1} 






{1,1,2} 




{4} 


{1,1,1,1} 






{1,1,2} - 


■> {1,3} - 


{4} 


{1,1,1,1} 


— )■ 


{1,1,2} 


^ {2,2} 


— > 


{4} 



(2.23) 



Note that parts of the map where the partition does not change (such as {1,1,2} — >• 
{1,1,2}) contribute nothing to the unitarity relation. 

So far we have only discussed phase space constraints from clustering. We will now 
include the effects of the measurement, where much of this structure is employed. For 
reference, we give a brief dictionary of symbols defined up to this point in Table 1. 



3. The Measurement Function for Clustering Algorithms 

The outcome of a jet algorithm is associated with a measurement on the final state. Each 
partition V gives a contribution A4('P) to the measurement function, where A^(P) is the 
Fourier transform of the measurement function. To find the total n-particle measurement 
function we can simply weight each term in the unitarity relation of Eq. (2.8) by the 
measurement factor for that partition, so that 



A^W = ^iV(P) J]c(p) 



V 



p&V 



M{V)s{V). (3.1) 



This is a more specific form of Eq. (1.7) in transform space. A similar form of this equation 
applies to momentum space; with the appropriate mapping Ai ^ M. The measurement 
function in Eq. (3.1) can be applied to any standard clustering algorithm with any observ- 
able. In order to simplify the effect of clustering on the structure of M^"^' we would like 
to use the maps developed in Sec. 2.2 to rewrite s(P) in terms of clustering and merging 
functions. In order to do this we require that the observables that we measure, {x} in 
Eq. (1.3), are only a function of momenta in each subset p and do not depend on multiple 
subsets in V. This is made explicit in Eq. (1.7), and hence in Eq. (3.1), where the phase 
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Symbol 


Description 


V 


Partition giving the outcome of the algorithm. 
See Eq. (2.4). 


c{p) 


Requires p particles to cluster. 
See Eq. (2.5). 


m{{pi,...,pk}) 


Requires groups oi pi,. . . ,pk particles to merge. 
See Eq. (2.13). 


s({pi,...,Pfc}) 


Requires groups oipi,...,pk particles to not cluster. 
See Eq. (2.6). 


-po 


The base partition {1, . . . , 1}. 
See Eq. (2.18). 


-pw 


The partition giving the outcome of the algorithm, used in maps. 

See Eq. (2.18). 


-pc 


A coarser partition than V^ , obtained by clustering single particles. 
See Eq. (2.23) for an example. 


rf 


The final step in a map, obtained from V^ by merging or clustering. 

See Eq. (2.18). 


MV'",v^) 


A map from V^ to V^, used in rewriting s{V^). 
SeeEq. (2.21). 


UV'.v^) 


A map from V^ to V-^ , used in rewriting s{V^). 
See Eq. (2.21). 



Table 1: Dictionary of symbols defined to describe phase space constraints from clustering. 

space constraints from clustering are separated from the measurement^. This means that 
the contribution to the measurement of a given outcome V of running the jet algorithm 
can be written as 

M{V) = n Hp) , (3.2) 



where 



peP 



w{p) : contribution of p clustered particles to the measurement. 



(3.3) 



Note that the single particle measurement function is M^^> = w{l), so in the absence of 
clustering 



M{{1,...,1}) 



A^fd) 



(3.4) 



As an example of a measurement that satisfies Eq. (3.2), consider running an inclusive 
jet algorithm and measuring the jet mass for jets with px > p^^^- Each p & V corresponds 
to a candidate jet produced by running the clustering algorithm, for which we measure its 



^In principle, one can define observables that can only be written as a function of all the momenta in an 
event, and cannot be disentangled into contributions from separate regions of phase space. However, such 
observables tend not to arise in practical applications. 
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mass mFj (or Xj in Fourier space). In this case Ai{V) becomes 



M{V)= JJexp 



-zx^0(pt(p)>p|"*) 5^9^ 



Qj 



»jep 



(3.5) 



where qi are the momenta of particles in a given partition p. Note that for this observable, 
the one particle measurement function is A^(l) = 1 for massless final states. 

For the very general class of observables satisfying Eq. (3.2), the n-particle measure- 
ment function in Eq. (3.1) becomes 



7W(") = ^Ar(p) l[c{p)wip) 



s{V), 



(3.6) 



where all partitions have n particles. We now show how the maps in Sec. 2.2 can be used 
to rewrite s{V) to express the full measurement function in a form that makes the effect 
of clustering explicit. 

We can think of the measurement as a map from the base partition V^ = {1, . . . , 1} 
to the outcome of the algorithm V^ : 



measurement map 4>m{'P j'P^) '■ 'P 



V' 



(3.7) 



Each map is associated with the clustering and measurement factor for the partition. The 
full phase space constraints includes a factor of s{V), but we define the contribution from 
the map (j)m without this factor. The s{V) factors will be rewritten using the maps (/>«, and 
(j)c from Sec. 2.2. For 3 particles, the measurement maps and contribution to Eq. (3.6) are 



</>r 



■pu 


— >• 


-pto 


contribution 


{1,1,1} 


— ;> 


{1,1,1} 


-5^(1)3 


{1,1,1} 


— ;> 


{3} 


c(3)w(3) 


{1,1,1} 


— ;> 


{1,2} 


3c{2)w{2)w{l) 



(3.8) 



The measurement map corresponds to Eq. (3.6) with V replaced by V^. The remaining 
work is to rewrite s{V'") in terms of inclusive constraints, which we use the clustering maps 
for. 

Putting the measurement map together with the clustering maps defined in Sec. 2.2, 
we have a complete map from the base partition to the final partition: 



ytot 



(po^ P/) ^ cPmir'',rn ■ 4>n,{v^,v') ■ M^", v^) 



(3.9) 



If we write all such maps, then we can determine the outcome of the algorithm. 

As an example, we give the measurement functions for the 2 and 3 particle cases. In 
Fig. 7 we use the diagrammatic notation of Sec. 2 to show the measurement maps V^ — )• V^ 
and clustering maps V^ -^ V^ ^ V^ for the 3 particle case. For n = 2, there are three 
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Vo 



measurement 
map 



V 



44 



w 



clustering 
map 



vf 



+ -^^ - -^^ 



■©H 



■^^ 



Figure 7: The maps to determine the measurement function for 3 particles are shown diagram- 
matically. The measurement map 7^° — )■ V^ corresponds to Eq. (3.8), while the clustering maps 
-pw _^ -pf YevfYite s{V^) using the unitarity relations displayed in Fig. 5. The contribution from 
each map is given in Eq. (3.12). 



maps. 



<^. 



tot 



{1,1} 
{1,1} 
{1,1} 



{1,1} 
{1,1} 

{2} 



{1,1} 
{2} 
{2} 



{1,1} 
{2} 
{2} 



contribution 

-c{2)w{lf 
c{2)w{2) 



(3.10) 



Therefore the total measurement function is 



7W(2)={^(l)2 + c(2)[{?(2)- 
= [Mil " 






(3.11) 



which defines AM^"^'. For n = 3, when the measurement maps from Eq. (3.8) are combined 
with clustering maps, which rewrite s{V^), there are seven maps. 



-'tot 



•po 
•pO 
•pO 

pO 
pO 
pO 
pO 



{1,1,1} 
{1,1,1} 

{1,2} 

{1,1,1} 

{1,2} 

{1,1,1} 

{3} 



{1,1,1} 
{1,2} 
{1,2} 
{1,2} 
{1,2} 

{3} 

{3} 



T'-' : contribution 

{1,1,1} : ^(1)3 

{1,2} : [w{in-3c{2)] 

{1,2} : 3c{2)w{2)w{l) 

{3} : Rl)3][-3c(2)][-m({l,2})] 

{3} : [3ci2)wi2)w{l)][-mi{l,2})] 

{3} : [wiin-ci3)] 

{3} : ci3)wi3) 

(3.12) 
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These combined to give the measurement function 

M^^) = w{lf + 3c(2) [^^(2) - w{lf]w{l) 

+ (c(3) [w{3) - Silf] - 3c(2)m({l, 2}) [w{2) - w{lf]w{iy^ 
= [M^^'>] ^ + 3A{(i) AA^(2) ^ ^^(3) _ (3_^3) 

This equation defines AA^*-^'. Note that for the 2 and 3 particle cases, a factor of w{k) 
for k > 1 comes with an accompanying —iJo{l)^. This is true in general, which we can 
show by making use of the separation between the maps (j)w and (j)c- This will lead to a 
simplification of the structure of the measurement function. 

Consider the set of maps 4>m{'P^,'P^)'4'wi'P^i'P'^) from 7^^ to V^ with a fixed T"^. These 
maps are generated by taking all possible V^ . For example, for n = 5 and V^ = {1,2,2}, 
these maps are (with V^ in bold for each map): 



■pO 


— >■ 


pc 


contribution 


{1,1,1,1,1} - 


■> {1,1,1,2} - 


■> {1,2,2} 


I5ci2)^w{lf 


{1,1,1,1,1} - 


■> {1,1,1,2} - 


■> {1,2,2} 


-30c{2)^wi2)wil)^ 


{1,1,1,1,1} - 


■> {1,1,1,2} - 


■> {1,2,2} 


15ci2fw{2fwil) 

(3.14) 



The sum of these terms is 



I5ci2f[wi2) - w{lf]wil) . 



(3.15) 



This binomial structure persists. All terms will share the same clustering factors, differing 
only in the combinatoric and measurement factors. If we sum over all V^ for a given V^ 
(V^lV^), then the clustering and measurement factors are 



'pruypc 



(3.16) 



p&v 



peV" 



The combinatoric factor is defined relative to N{V'^), and so it is the number of ways to 
distinguish the map V^ — )• V^ — )• V'^ from V^ — )• V'^. This is simply 



NUV'",V') 



m 

k=2 ^^k 



(3.17) 



where d^ is the number of elements of V^ equal to k (recall Eq. (2.3)). Therefore, we can 
sum over all V^ and Eq. (3.16) becomes 



NiV) \ Yl c{p)Aw{p) 



where 



Aw{p) 



w{p) - iB{l)P 
w{l) 



if p > 1 , 
if p = 1 . 



(3.18) 



(3.19) 
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Finally, we must deal with the map (j)ciJ-'^,'P^). Given a partition V^, we will define 
the set of all maps (pc{'P'^,'P^) (all such maps for all allowed V^) as Qc{V^). For example, 

g,{{l, 1, 2}) = {(/.e({l, 1, 2}, {1, 1, 2}), (/>e({l, 1, 2}, {1, 3}), ,/.e({l, 1, 2}, {4})} . (3.20) 

These maps from V^ to V^ are the second, third, and eighth maps in Eq. (2.23). For each 
map (pc, we define 

J-{(f>c) ■ the product of inclusive constraints for the map (pc ■ (3.21) 

Using the maps in Qci{l, 1, 2}) as an example, 

^((/<e({l,l,2},{l,l,2})) = l 

J-((/),({l,l,2},{l,3})) = -2m({l,2}) 

T{M{h 1, 2}, {4})) = 2 m({l, 2}) m({l, 3}) . (3.22) 

3.1 Main Results 

Putting the pieces together, the measurement function is 






■pc 



pev 



(3.23) 



This is one of our main results, and can be used to arrange the measurement function into 
a form suitable to study the logarithmic contributions to the cross section. 

Each map (j) : V^ — )■ V^ will contain some number of single particles not participating 
in any clustering or merging. These particles will only be associated with factors of w{\) = 
M^^'. Since M.^^' only depends on the phase space constraints from a single particle, it is 
useful to distinguish these contributions from multi-particle phase space constraints^. 

Suppose we consider maps with exactly k factors of M^^'. These single particles do 
not participate in the maps at all, and will remain in the final partition V' for each map. 
Therefore we can factor out these particles from the map, and they will become maps on 
n — k particles: 

(pn -^ 4>n-k (3.24) 

If the map (pn leads to a contribution 7W ((/)„) to the measurement function, then we can 
write this as 



M{(Pn) 



n\ 



-Mi 



^n—k. 



[M^ 



1)1 fc 



(3.25) 



{n-k)\k\ 

The combinatoric factor changes the 1/n! from the n-particle matrix element into separate 
factors for M{(t)n-k) and [A^^^^] . This gives another of our main results. 



n-2 



Al(«)=[A^(i)]" + ^ 



n\ 



fc=o 



{n-k)\k\ 



A7W("^^)[A^ 



■{1)1 fc 



(3.26) 



*When we consider boundary clustering, in many cases M^^' is the only factor that contributes at 
leading log to the cross section. For an example, see [27], where the contribution AtM"^ was calculated 
for the C/A and kx algorithms for a jet mass observable in dijet events. This will be discussed further in 
Sec. 5. 
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where AA^*^"^ are the n-particle measurement function terms where all particles cluster or 
merge with at least one other particle, and we can use this equation to define ATW'"-*. For 
2 particles, 

AM^^^ = c{2)[w{2) - iBilf] 

= c{2)Aw{2) , (3.27) 

and for 3 particles, 

A7M(3) = c(3) [w{3) - w{lf] - 3 c(2)m({l, 2}) [w{2) - w{lf] w{l) 

= c{3)Aw{3) - 3 c(2)m({l, 2})Aw{2)w{l) . (3.28) 

We give the measurement functions for n = 4 and n = 5 in Appendix A as additional 
examples. 

4. Cancellation of Infrared Divergences from Clustering 

The formalism for the measurement function lets us study the IR structure of the various 
clustering contributions for the Abelian terms. As an example, we consider the 3 particle 
terms in the measurement function in the limit that one particle is very soft. This argument 
can be generalized to higher order terms, or adapted to the collinear limit. For 3 particles, 

A^(3) = [M^Y + 3A>('(2)_^(i) ^ ^^{-i) 

= w{lf + 3c(2) [w{2) - w{lf]w{l) 

+ c(3) [w{3) - w{lf] - 3c(2)m({l, 2}) [w{2) - w{lf]w{l) . (4.1) 

The single particle measurement function A4^^' leads to an IR safe contribution to the cross 
section. We will consider AA^^^\ which is the last line of Eq. (4.1); in showing these terms 
are IR safe, we will show that terms involving AA^'^^ are also free of IR divergences. 

Let us label the soft momenta ki, k2, and ^3; we will take the limit ^3 — )• 0. The term 
involving the constraints c(2)m({l, 2}) has been symmetrized over soft momenta, meaning 
if we want to take ^3 soft then we must split this up into 3 separate terms. Labeling the 
momenta, this term becomes 

3c{2)m{{l,2})[w{2) - w{lf]w{l) 

-^ c{{ki, k2})m{{k3, {^1,^:2}}) [w{ki, ^2) - w{ki)w{k2)]iij{k3) 
+ c{{ki, k3})m{{k2, {ki, k-s}}) [wiki, ks) - w{ki)w{k3)] w{k2) 
+ c{{k2, k3})m{{ki, {k2, ks}}) [w{k2, ks) - w{k2)w{k3)] w{ki) . (4.2) 

In this language, the IR safety of AA^'^^ is due to the relations 

lim w(k, kg) = w{k) , lim w{ks) = 1 . (4.3) 

This implies that 2 particle clustering effects are IR safe: 

lim w{k, ks) - w{k)w{ks) = => lim AA^^^) ^ g . (4.4) 
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This relation simplifies the terms in Eq. (4.2) to 

c{{ki, k2})m{{k3, {ki, /C2}}) [w{ki, k2) - w{ki)w{k2)\ ■ (4.5) 

The relation in Eq. (4.3) can be generalized to 

lim w{n) — )• w{n — 1) , (4.6) 

which means that the c(3) term in ATW^"^) becomes 

C{{ki,k2, ks}) [w{{ki, k2}) - W{ki)w{k2)] . (4.7) 

Therefore both groups of terms are IR divergent in the soft limit. However, by considering 
the merging constraints one can show for an IR safe jet algorithm that 

lim c({A:i, A;2})m({A;3, {ki,k2}}) = c{{ki, k2, hs}) . (4.8) 

The reason is that when ^3 is soft, clustering with any other particle ki would leave the 
momentum of ki unchanged. In this limit the constraints requiring the soft particle to merge 
with a clustered pair is the same as all three particles clustering. Therefore, both terms in 
AA^(^) reduce to Eq. (4.7), and they cancel (AA^*-"^' — )■ 0). This is an explicit demonstration 
that IR divergences are not introduced by clustering effects for 2 or 3 particles. 

5. An Application: All-Orders Abelian Structure of Soft Clustering for 
Jet Shapes 

The all-orders description of the measurement function developed in the previous sections is 
useful because it enables us to understand the effects of clustering on the higher order terms 
in the cross section for jet observables. As an example we consider e'^e~ — )• N jets defined 
by an inclusive jet algorithm (such as C/A) and define two observables, p and A, where 
p = (ml + . . . +'m'j^)/Q^ is the sum of jet masses {rrii} scaled by the center of mass energy 
Q and A is the total energy of radiation outside of all jets. When {Qp,A} <^ Q^/p ^ Q, 
each jet mass ttt-j is forced to be small and the event is described by collinear radiation 
in the jets and soft radiation in and between the jets. This observable was considered for 
dijet events in [30, 26, 27, 31, 32, 33]. The all-orders structure of the Abelian terms in this 
section apply to a wide range of jet shapes in the soft/collinear regime, such as angularities. 
Perturbative calculations generate logarithms of the form Inp and InA/Q, and a reli- 
able prediction in this regime requires resummation. This is achieved through a factoriza- 
tion theorem, for which for the double cumulant takes the form 



where 



j2 

aoHNiQ) Mp) ^ ■ ■ ■ ^ Jn{p) ®57v(p,A). (5.2) 



dpdA 
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The hard function H is determined by the short distance interaction that produces the high 
energy jets, and depends on the center of mass energy Q. The jet functions parameterize 
cohinear radiation for each jet at the scale Q^/p and the soft function parameterizes soft 
radiation in and between jets at the scale Qp. The radiation collinear to each jet direction 
clusters only within the jet when y^ <^ R, and in this case the effect of clustering is 
described entirely by the boundary clustering of soft particles in the event, as discussed in 
Sec. 1. It was shown in [26, 27] that boundary clustering gives rise to clustering NGLs in 
the soft function of the form ln{Qp/A) starting at 0{a'^ In ). 

Sn{p,A) = ^{0\WlyMi{p,A})W{ny\0) (5.3) 

where W/^} represents a product of soft Wilson lines, Yj^, where n^ = (1, n,) is a light like 
vector in the i-th jet direction, with the appropriate color representation and path ordering 
of the final state hard partons^. The soft function is a matrix in color space and cat is the 
color normalization factor such that S]\f{p,A) = 5{p)6{A) at tree-level. The measurement 
operator Ai {p, A) implements the measurement of p and A for a given jet algorithm, and, 
for final states \Xs) with soft momenta {k}, the measurement function acts as 

Mip,A)\Xs)=M{p,A,{k})\Xs). (5.4) 

For convenience we work with the Fourier space soft function, 

/oo 
dpdA exp{—ixp) ex-p{—iyA/Q)S]\[{p, A) , (5.5) 

-oo 

for which the measurement function is typically expressed in the form 

M{x,y,{k}) =Y[expl-ixY,'^^T^J.{kj-Ak}) exp -iy ^ ^7^out(A:,•; {/c}) 

(5.6) 
where the functions TZj^{kj\ {/c}) and 7^out(^j; {^}) implement the phase space constraints 
of the jet algorithm, requiring particle kj to be in the i-th jet and out-of-jet regions re- 
spectively. Each region, TZi, depends on all the soft momenta in the final state {fe} and 
therefore implicitly on all of the other regions. The formalism developed in this paper helps 
disentangle this dependence. For n final state particles it is useful to organize the effects 
of clustering as in Eq. (1.13). When the measurement function is integrated against the 
squared matrix element from the Wilson lines, we get the ©(a") contribution to the soft 
function, 

SW({X}) = j fn^) -4W({A:})A?W({x},{A.}) . (5.7) 

^Recall a soft Wilson line directed along ri'', in representation T", contains soft gauge fields (AJ^ ) and is 
defined by a path ordered exponential: 



Y^{x) = Pexp [ig I ds n ■ A'^{x + sn) T' 
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The first term in Eq. (1.13) is the measurement function for an algorithm for which the 
measurement function factorizes into single particle measurements in the soft sector, such 
as for the anti-kx algorithm. The second term accounts for correlations in the measurement 
function. This leads us to write 

5W({x}) = S^^{{x}) + Si±{{x}) . (5.8) 

The algorithm-dependent effects of soft clustering are entirely contained in the last term. 
When soft gluons boundary cluster they can not be collinear to the jet direction, 
as depicted in Fig. 3. In general a final state gluon can contribute a double log to the 
cumulant cross section, ag In (Qpc/Ac), associated with soft and collinear divergences. For 
the Abelian terms in the soft matrix element, collinear divergences arise only for soft gluons 
collinear to the jet direction. Therefore, as shown in [27], those gluons involved in boundary 
clustering can each contribute at most a single log, Os ln((5/9c/Ac), from a soft divergence 
to the Abelian terms in 5corr. For example, at 0{a^), A4cott contains a contribution where 

~(3) 

all three particles boundary cluster together, which contribute at NLL to 5'corr. However 
the measurement function also contains the contribution 

^A^' ' D 6 (gluon 1 and 2 boundary cluster) 

X 9 (gluon 3 is in the jet and does not boundary cluster) , (5.9) 

which a priori could contribute at LL in the exponent, a" In"* for m > n, since gluon 3 is 
sensitive to both soft and collinear divergences. In order to show that clustering NGLs do 
in fact occur at NLL, we must show that constraints of the form Eq. (5.9) can be expressed 
in terms of lower order phase space constraints. Rewriting the measurement function in 
terms of inclusive constraints allows us to do precisely this. For the term in Eq. (5.9) this 
is illustrated graphically is Fig. 2. Next we will use our main result, Eq. (3.26), to show 
the all-orders form of the Abelian soft function. We will see that clustering NGLs give new 
contributions at NLL at all orders, spoiling the exponentiation of the Abelian terms. 



5.1 All-Orders Abelian Structure 

A considerable simplification happens in the Abelian sector. It is well known that the 
(D{a^) Abelian matrix element factorizes as 

1 " 

A^"\{ks}) = -J[AW{k,), (5.10) 

where A^^'{k) is the squared matrix element at 0{as). When clustering effects are absent 
or power suppressed, the measurement function factorizes into the product of single gluon 
constraints, as for the anti-kx algorithm in Eq. (1.11). In this case it is easy to see that 
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Abelian soft function exponentiates. The 0{a^Cp) contribution to the soft function is' 

1 / /■ d^A; 



'Abel 



(W) 



1 



(2^)' 



A^^\k) M^^\{x}) 



5«({x}) 



which is summed to produces the-ah orders Abehan soft function, 



SAbei.iix}) = exp 



5«(W) 



(5.11) 



(5.12) 



This is a manifestation of Abehan exponentiation and we see that factorization of the 
measurement function, as in Eq. (1.11), is a necessary condition for it to hold. 

For algorithms for which the measurement function contains corrections from cluster- 
ing, we can use our main result, Eq. (3.26), to express these corrections as 



Mtl{{^},^)=Y. 



n\ 



k=0 



{n-ky.k\ 



A7W ("-'=) [A^(i 



(5.13) 



where AA^'*^-* is the part of the n-particle measurement function in which all n particles 
cluster with at least one other particle. These corrections modify Eq. (5.11) to give 

n-2 



^il. = ;^(^^^^)" + E^^£-'^^(^^^^)', 



(5.14) 



fc=0 



where we have defined 



A5irJ(W) = i,ny |^^('Hfc.)AA^W({x},{A:}). (5.15) 

Finally, although the Abelian soft function no longer exponentiates, the form of Eq. (5.14) 
allows us to write the all orders soft function as 



5Abei.({x}) = exp 5«({x}) l + J^A^lgdx}) 



(5.16) 



k=2 



5(fc) 



Since AS'g^J arises from phase space constraints where each of the k gluons boundary 
clusters with at least one other gluon, it contributes at most at NLL; that is aj In' {Qp^/ A) 
in the cumulant cross section. Resummation of clustering NGLs would require a relation 

~(n) 

between AS'g^j for all n. However, as we have shown from the general form the measurement 
function in Eq. (3.23) and Eq. (3.26), at each order in a" the clustering terms in the soft 
function receive a contribution from 



AA^(") D c{n)[win) - wiiy""] 



(5.17) 



where all n particles cluster with each other. This contribution is new and unique at 
this order in a^. This makes resummation of clustering NGLs very unlikely. In the next 
subsection, we will see that the all-orders form of the Abelian soft function in Eq. (5.16) is 
required in order to have a renormalizable soft function. 

®In pure dimensional regularization, Abelian matrix elements with a virtual gluon are scaleless, hence 
we only consider the matrix element where every particle is in the final state. 
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5.2 Renormalization of the Abelian Soft Function 

It is well known that the renormalization of the product of Wilson lines, such as those that 
appear in the soft function, obey Abelian exponentiation - that is, the Abelian contribution 
to the all-orders anomalous dimension depends only on the the one-loop result. However, 
the presence of the measurement operator in the soft function makes its structure and 
renormalization more complex. We now consider the effects that clustering has on the all 
orders UV structure of the Abelian soft function. 

We have seen in Eq. (5.12) that since the anti-kx measurement function can be ex- 
pressed as the product of single gluon constraints, the soft function in this case obeys 
Abelian exponentiation. This means that the bare anti-kx soft function is related to the 
renormalized one to all orders as follows: 






-^akT "S^akT- 



exp 



?(1) 



(5.18) 



where b stands for bare, r stands for renormalized, and Z^^ is the Abelian position space 
one- loop counter-term. We have suppressed the subscript Abel.' since all of our statements 
in this section only refer to Abelian terms. Consistency of the factorization in Eq. (5.2) and 
the requirement that the cross section must be renormalization scale independent imply 
that the anomalous dimensions must satisfy 



N 



IH + Yl 7.h + IS 



0, 



(5.19) 



Since the hard and jet function anomalous dimensions do not depend on the algorithm, 
the soft function anomalous dimension must also be independent of it. This implies that 



?(i) 



?(i) 



Zg^j = Z^-^ must hold for all algorithms. Any contribution to the soft function AS'aig must 
be UV finite, since no new divergences are introduced by these algorithms when compared 
to anti-kx. 

Since algorithms in the kx class share the same Z factor, the bare soft function is 
related to the renormalized one as follows: 



'-'alg ~ ^akT "^alg ~ ^^P 



z 



(1) 

akT 



qr 
'-'alg ' 



(5.20) 



using Eq. (5.18). When S^^' is replaced with its bare version in Eq. (5.14), it takes the 
form ^ 



;6(n) 



1 



n\ 



(5^1))" + ^A5(' 



n—k) 



k=0 



k\ 



(5.21) 



We can rewrite this to all orders as 



5.^ 



alg 



exp 



5^1)1 1 + ^A5, 



(fc) 

alg. 



(5.22) 



fe=2 



After inserting S^^^' = Z^^ + S^'^^> into the previous equation, it follows from Eq. (5.20), 
that the renormalized soft function is 



5:ig = exp[s'^WJ 1 + X^A5, 



(5.23) 



fc=2 



27 



~(n) 

Lastly, we also know from Sec. 4 that AS'^ is IR finite in addition to UV finite. As 
such, these terms can be calculated numerically without the need for difficult subtraction 
schemes. This is a major advantage, since the phase space constraints imposed by AA^'") 
are typically quite complex, making the effects difficult to calculate analytically. The O(a^) 
corrections for the C/A and kx algorithms were calculated in [27] for the dijet case {N = 2). 

5.3 Failure of Exponentiation 

We might hope that the corrections due to clustering in the measurement function could be 
redefined to reveal further structure in the measurement function beyond that in Eq. (3.26). 
In particular, one might expect that the corrections to the measurement function could be 
redefined in such a way that the Abelian corrections exponentiate: 



•S'alg = exp 



5'r-(i)_^^^5"W 



(5.24) 



=2 

This is equivalent to reorganizing the measurement function to define a new correction 
term, ATW'^"' that is integrated against the Abelian matrix element to get AS"("\ as in 
Eq. (5.15). The form in Eq. (5.24) follows from the properties of clustering in the measure- 
ment function if we can define AS"'""^ at each order to involve only the clustering/merging 
of all n particles with each other. This is contrast with the form in Eq. (5.23), where AS^^' 
involves clustering effects of n or fewer particles. Of course, we could simply define AS"' 
through Eq. (5.24) at each order, since the correction terms AS^^J in Eq. (5.23) are finite. 
However, unless we show Eq. (5.24) from the properties of clustering in the measurement 
function, it has no more content than Eq. (5.23). 

Terms grouping all n particles together arise, in the language of Sec. 3, from mappings 
that have V^ = {n}. To illustrate the difference with an explicit example, consider n = 4, 
which is worked out in full detail in Appendix A. AS^'^' contains clustering effects of 4 
particles as well as clustering effects from 2 particles in two separate groups: 

AM^^^ = 3 [A>f(2)]2 ^ ^_y^' (4) ^ (5_25) 

In Appendix A, we see that AA4^^^ can be rewritten such that every term in AA^'^^^ 
merges or clusters all 4 particles into one group and is given in Eq. (A. 7). In general, this 
approach requires defining AA^ ' ("^ such that clustering effects at lower orders (ATW ' ^''' 
for k < n) are explicitly factored out. However, for n = 5 (and higher), this is not the case, 
as demonstrated in Appendix A. With 5 particles, there is a leftover term that cannot 
be placed into products of AM '^^' for fc < 5, nor does it cluster or merge all 5 particles 
together. This term must reside in AM'^^>: 

AA{'(5) :^ 2,{)c{2fm{{l,2})[w{2) - w{lf]w{lf . (5.26) 

It contains two groups, one of 2 particles and one of 3 particles. Such terms are associated 
with the map from V'^ = {1, 1, 1, 2} — t- P' = {2, 3}, when clustering takes place after the 
first merging step. For n > 5, there are terms contributing to AA^''"' in which each of 
the n particles cluster or merge with at least one other particle but all n are not grouped 
together, i.e. V^ ^ {n}. We are forced to conclude that Eq. (5.24) does not have any 
content and is simply a redefinition of Eq. (5.23). 
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6. Conclusions 

In this paper, we developed a formalism to determine the phase space constraints from 
a clustering algorithm and the effects on a measurement function. This framework can 
be used to study perturbative corrections that arise from clustering. We expressed the 
constraints from clustering as a function of the outcome of the algorithm. This is charac- 
terized in terms of a partition that specifies which particles are clustered by a clustering 
function, c{p), and a non-clustering function, s{{p}), which requires groups of particles not 
to cluster. We derived a set of unitarity relations on phase space, which can be iteratively 
applied to eliminate the non-clustering function. As discussed in the text, this requires the 
introduction of a merging function, m{{p}), which appears in unitarity relations involving s 
functions. This rewriting eliminates exclusive phase space constraints that require particles 
not to be in certain regions of phase space, producing only inclusive constraints. This is 
useful because it simplifies the structure of the most divergent regions of phase space and is 
a key step in allowing us to relate higher order clustering effects to lower order constraints. 

In Sec. 3, we applied this result to the measurement function, weighting the phase 
space constraints for each outcome of the algorithm by the contribution to the measure- 
ment. Grouping terms that have the same phase space constraints from clustering, we 
rewrote the measurement function purely in terms of inclusive constraints in Eq. (3.23), 
eliminating the non-clustering function s. This form makes the clustering effects explicit, 
expressing the measurement function as a term associated with no clustering and a set of 
corrections. These correction terms are associated with n-particle clustering effects, and 
are proportional to factors defined in Eq. (3.19) that are the difference between the contri- 
bution to M with and without clustering. This made the analysis in Sec. 4 of the infrared 
structure associated with clustering more straightforward. Additionally, we showed that 
the n particle measurement function A^(") can be related to the measurement function 
appearing at lower orders, Ai^^' with A; < n, as in Eq. (3.26). This form is particularly 
useful when trying to understand the higher order structure of jet cross sections. 

In Sec. 5 we apply our formalism for the measurement function to examine the per- 
turbative corrections that arise from clustering for a wide class of jet shapes in the soft 
collinear regime. The main results for the measurement function, chiefly Eq. (3.23) and 
Eq. (3.26), were used to determine the all-orders structure of the Abelian terms in the soft 
function, Eq. (5.23). This structure is universal to jet shape observables and depends only 
on the form of the measurement function and constraints from renormalization. From the 
all-orders form of the soft function and by count logarithms arising from clustering in the 
Abelian terms, we found that to all orders clustering NGLs contribute at NLL plus sub- 
leading logs for the kx class of algorithms. This same framework can be applied to study 
non- Abelian terms or non-logarithmic contributions from clustering, since the measurement 
function contributions from clustering are made explicit. 
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A. Examples for 4 and 5 Final State Particles 



In this appendix we present the measurement function for 4 and 5 final state particles using 
the formalism of Sec. 2 and Sec. 3. We will start with 4 final state particles, for which 
there are five possible partitions. 



{{4}, {1,3}, {2, 2}, {1,1, 2}, {1,1, 1,1}}. (A.l) 



As dictated in Eq. (3.23), we split the contribution into two factors: the one associated 
with V^ and the one associated with the merging sequence. For each allowed V^, the 
appropriate factors from the first half of Eq. (3.23) are: 



P^ 


: N{V'')Y[^.^T,.c{p-)/^w{p-) 


{4} 


: c{4:)Aw{A) 


{1,3} 


: 4:c{3)w{l)Aw{3) 


{2,2} 


: 3c(2)2A-55(2)2 


{1,1,2} 


: 6c{2)w{lfAw{2) 


{1,1,1,1} 


w{l)^. 



(A.2) 



Each of these factors is then multiplied by X](Aeg„(pc) ^{4>)- The function T{(p), defined in 
Eq. (3.21), depends on the specific sequence of mergings that is represented by the mapping 



: V' = Vf ^Vt--- rf -^ rf^i = Vf . (A.3) 



Recall, from the definition of this map, that the first step in this mapping, Vf — s- V2, must 
be a merging step and not a clustering clustering. All of the possible mappings for n = 4 
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and the associated factors are listed below. 



vf 


H^) 


{4} 


1 


{4} 


-m({l,3}) 


{4} 


-m({2,2}) 


{4} 


-m({l,l,2}) 


{4} 


2m({l,2})m({l,3}) 


{1,3} 


1 


{1,3} 


-2m({l,2}) 


{2,2} 


1 


{1,1,2} 


1 


1,1,1,1} 


1 



ct>a : P" > 

{4} > 

{1,3} > 

{2,2} > 

{1,1,2} > 

{1,1,2} -^ {1,3} -^ {4} : 2m({l,2})m({l,3}) (A.4) 

{1,3} > 

{1,1,2} > 

{2,2} > 

{1,1,2} > 

{1,1,1,1} ^ 

After combining the factors from Eq. (A.4) with appropriate factor from Eq. (A. 2), as 
dictated by Eq. (3.23), we get 

A^(^) = w{lf + c(4)AiI^(4) + 4c(3)5J(l)At2(3)[l - m({l, 3})] 

+ 3c(2)2Aw(2)2[1 - m({2, 2})] + Qc{2)w{lf /:^w{2)[l - m({l, 1, 2})] 

- I2c{2)w{lfl\w{2)m{{l, 2})[1 - m({l, 3})] (A.5) 

We can write the expression for 7W'^^ in the form of Eq. (3.26) by substituting the expres- 
sions for AX (2) and AA^^^) given in Eq. (3.27) and Eq. (3.28). The resuh is 

where 

AA^(^) =c(4)Aw(4) - 4c(3)iy(l)Atiy(3)m({l, 3}) + ?,c{2f^w{2f \l - m({2, 2}) 

+ 6c(2)tI5(l)2A^5(2) 2m({l,2})m({l,3}) -m({l,l,2}) . (A.7) 

Next, we give the primed version as defined under Eq. (5.24): 

AaT'(^) = AAT^^) - 3[A>('(2)]2 (A_8) 

= c(4)Aw(4) - 4c(3)w(l)AzI5(3)m({l, 3}) - 'ic{2fh.w{2fm{{2, 2}) 

+ 6c(2)'5;(l)2AiI^(2) 2m({l,2})m({l,3})-m({l,l,2}) . (A.9) 

Each term is associated with a map with V^ — )• {4}. 
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For n = 5, we proceed similarly. There are seven possible partitions, each contributing 
to the measurement function: 

{{5}, {1,4}, {2,3}, {1,1, 3}, {1,2, 2}, {1,1, 1,2}, {1,1, 1,1,1}}. (A.IO) 

The factors associated with clustering are listed below. 



pc 


: N{V^)Up'^^r^c{p-)Aw{p-) 


{5} 


: c(5)Aw(5) 


{1,4} 


: 5c{4:)w{l)Aw{4) 


{2,3} 


: 10c{3)c{2)Aw{2)Aw{3) 


{1,1,3} 


: l0c{3)w{l)^Aw{3) 


{1,2,2} 


: 15ci2)^wil)Aw{2f 


{1,1,1,2} 


: 10c{2)w{lfAw{2) 


{1,1,1,1,1} 


wil)^ 



(A.U) 



Next we give the factors associated with merging/clustering. For ease of presentation, 
we give those maps with V-' = {5}. 



(piT'' 


> 




p/ , 


n<t^) 


{5} 


>■ 




{5} : 


1 


{1,4} 


> 




{5} : 


-m({l,4}) 


{2,3} 


> 




{5} : 


-m({2,3}) 


{1,1,3} 


> 




{5} : 


-m({l,l,3}) 


{1,1,3} 


-^ {1,4} 


— > 


{5} : 


2m({l,3})m({l,4}) 


{1,2,2} 


> 




{5} : 


-m({l,2,2}) 


{1,2,2} 


-^ {2,3} 


— > 


{5} : 


2m({l,2}m({2,3})) 


{1,2,2} 


-^ {1,4} 


— > 


{5} : 


m({2,2})m({l,4}) 


{1,1,1,2} 


> 




{5} : 


-m({l,l,l,2}) 


{1,1,1,2} 


-^ {1,4} 


— > 


{5} : 


3m({l,l,2})m({l,4}) 


{1,1,1,2} 


-^ {1,1,3} 


— > 


{5} : 


3m({l,2})m({l,l,3}) 


{1,1,1,2} 


-^{1,1,3}^{1 


,4}^ 


{5} : 


-6m({l,2})m({l,3})m({l,4}) 


{1,1,1,2} 


-^{1,1, 3} ^{2, 3}^ 


{5} : 


-3c(2)m({l,2})m({2,3}) 



(A.12) 
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The rest of the maps are now given. 

0c : P" > 

{1,4} > 

{1,1,3} > 

{1,2,2} > 

{1,1,1,2} > 

{1,1,1,2} - 
{2,3} 

{1,2,2} 

{1,1,1,2} - 

{1,1,3} 

{1,1,1,2} 

{1,2,2} 

{1,1,1,2} 

{1,1,1,1,1} 



{1,3} 

> 

> 

{1,1,3} 

> 

> 



vf 


^{^) 


{1,4} 


1 


{1,4} 


-2m({l,3}) 


{1,4} 


-m({2,2}) 


{1,4} 


-3m({l,l,2}) 


{1,4} 


6m({l,2})m({l,3}) 


{2,3} 


1 


{2,3} 


-2m({l,2}) 


{2,3} 


3c(2)m({l,2}) 


{1,1,3} 


1 


{1,1,3} 


-3m({l,2}) 


{1,2,2} 


1 


{1,1,1,2} 


1 


{1,1,1,1,1} 


1 



(A.13) 



We can construct M^^' by combining these factors using Eq. (3.23). After using the 
definitions of AA?^^)^ /\M^^) , and AA4(^), we get 

M^^) = [A^(i)]5 + 10 A>('(2)[_^(i)]3 ^ 10 AAT(3)[_y^{i)]2 + 5 A>((4)_y^(i) ^ ^_^(5) ^ 

(A.14) 

where 

A7W(5) = c(5)AiIJ(5) - 5c(4)tIJ(l)Aw(4)m({l,4}) 
+ 10c(3)c(2)AiZ?(2)A{iJ(3) l-m({2,3}) 
+ 10c(3)ty(l)2Aziy(3) [2m({l, 3})?n({l, 4}) - ?n({l, 1, 3}) 
+ lf>c{2fw{l)I\iB{2f m({l, 4})m({2, 2}) - 2m({l, 2}) (l - m({2, 3})) - m({l, 2, 2}) 
+ 10c(2){?(l)^A^I5(2) r-3c(2)m({l, 2})m({2, 3}) + 3m({l, 4})m({l, 1, 2}) 

+ 3m({l, 2}) (-27n({l, 3})7n({l, 4}) + m({l, 1, 3})) - ?n({l, 1, 1, 2}) 
+ 2,{)c{2fw{lf /\w{2)m{{l, 2}) . (A.15) 
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The primed version is: 

= c{5)Aw{5) - 5c{i)w{l)Awi4:)m{{l,4}) - 10c(3)c(2)A{5(2)Aw(3)m({2, 3}) 
+ Wc{3)w{lfAw{3) 2m({l, 3})m({l, 4}) - m({l, 1, 3}) 

+ l5c{2fw{l)Aw{2f [m({l, 4})m({2, 2}) + 2m({l, 2})m({2, 3}) - m({l, 2, 2}) 
+ Wc{2)wilfAw{2) -3c(2)m({l, 2})m({2, 3}) + 3m({l, 4})m({l, 1, 2}) 
+ 3m({l, 2}) (-2m({l, 3})m({l, 4}) + m({l, 1, 3})) - m({l, 1, 1, 2}) 
+ 30c(2)^io(l)^A'ji;(2)m({l, 2}) . (A.16) 

All of the terms in ATW^^-* are associated with maps that have V-^ = {5}, except for 
the last one. This term, whose map is shown in bold in Eq. (A. 13), is associated with a 
map from V^ = {1, 1, 1, 2} to V-' = {2, 3}. This contribution is associated with two groups 
of gluons, one with 2 gluons and the other with 3, rather than from clustering of all 5 
gluons. This is the first order at which clustering can occur after a merging step, and it is 
these terms which cause the difficulty. 
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